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Glioblastoma remains one of the most lethal types of cancer, and is the most 
common brain tumour in adults. In particular, tumour recurrence after surgical 
resection and radiation invariably occurs regardless of aggressive chemotherapy. 
Here, we provide evidence that the transcription factor ZEBl (zinc finger E-box 
binding homeobox 1) exerts simultaneous influence over invasion, chemo- 
resistance and tumourigenesis in glioblastoma. ZEBl is preferentially expressed 
in invasive glioblastoma cells, where the ZEBl-miR-200 feedback loop 
interconnects these processes through the downstream effectors ROBOl, c- 
MYB and MGMT. Moreover, ZEBl expression in glioblastoma patients is 
predictive of shorter survival and poor Temozolomide response. Our findings 
indicate that this regulator of epithelial-mesenchymal transition orchestrates 
key features of cancer stem cells in malignant glioma and identify ROBOl, 0LIG2, 
CD133 and MGMT as novel targets of the ZEBl pathway. Thus, ZEBl is an 
important candidate molecule for glioblastoma recurrence, a marker of invasive 
tumour cells and a potential therapeutic target, along with its downstream 
effectors. 



INTRODUCTION 

With a median survival of about 15 months [Schwartzbaum 
et al, 2006; Stupp et al, 2005), glioblastoma is the most frequent 
and aggressive of all gliomas, with a propensity to invade the 
surrounding parenchyma [Wen & Kesari, 2008). Individual 
tumour cells can be found far from the primary tumour site, 
often crossing great distances into the contralateral hemisphere 
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[Wilson, 1992). These cells cannot be isolated for surgical 
resection, or easily targeted by irradiation, and thus represent 
sources for tumour recurrences [Glas et al, 2010). Adjuvant 
chemotherapy [e.g. Temozolomide, TMZ) is therefore included 
as a critical component of the current standard of care, in attempt 
to address these residual invasive cells. Given the exceedingly 
poor prognosis, it is critical to understand the biology of 
treatment-resistant glioblastoma cells. 

The cancer stem cell hypothesis [Reya et al, 2001) postulates 
an intra-tumoural hierarchy, where a small population of 
tumour cells has greater abilities to initiate and propagate 
tumours [Ignatova et al, 2002), rendering cancer stem cells an 
important therapeutic target [Vescovi et al, 2006). Cancer stem 
cells have been shown to be more invasive and therapy resistant 
than other cells of the same tumours [Bao et al, 2006; Cheng 
et al, 2011; Lathia et al, 2011). Tumour heterogeneity is a direct 
implication of the cancer stem cell hypothesis, and indicates 
that cell populations with different properties [such as drug 
resistance or higher capacity for tumour/recurrence formation) 
exist within the same tumour [Siebzehnrubl et al, 2011). The 
triad of tumourigenesis [cancer stemness), invasion and therapy 
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resistance is a fatal combination if merged in a single cell 
population, and renders such a population an important 
contributor to poor outcome. 

This triad is induced by Epithelial-Mesenchymal Transition 
(EMT) in cancers outside the CNS, where EMT is the major cause 
of invasion and metastasis (Chaffer & Weinberg, 2011), and 
cancer cells undergoing EMT have been shown to acquire stem 
cell traits and are frequently more therapy resistant (Mani 
et al, 2008; Polyak & Weinberg, 2009; Singh & Settleman, 2010). 
Therefore, EMT can generate cell populations that combine 
these three above-mentioned hallmarks. However, the role of 
EMT and related processes in brain cancer has received little 
attention thus far (Kahlert et al, 2012; Lu et al, 2012; Mikheeva 
et al, 2010), likely because the brain is lacking critical tissue 
components {i.e. epithelium and mesenchyme). Yet, it is 
conceivable that key invasion pathways overlap between CNS 
and other cancers, and that factors inducing EMT outside the 
brain also activate the triad of invasion, stemness and chemo- 
resistance in malignant gliomas. 

Many factors, including reduced cell adhesion [Asano 
et al, 2004), reduced matrix adhesion (Nakada et al, 2007), 
matrix protease secretion (Rao, 2003) and cytoskeletal remodel- 
ing [Giese et al, 1996) have been advanced as determinants of 
glioma invasion; several of these pathways are induced by EMT 
outside the CNS (Chaffer & Weinberg, 201 1) . Therapy resistance 
in glioma is mediated by expression of DNA repair enzymes (Bao 
et al, 2006), and/or expression of drug efflux transporters 
(Bleau et al, 2009) . Of particular note, 0-6-Methylguanine DNA 
Methyltransferase (MGMT) confers resistance to the standard 
of care drug TMZ (Bocangel et al, 2002). While MGMT is also 
expressed in several non-CNS cancers (Gerson, 2004), it is 
currently unknown whether EMT can induce MGMT expression. 
Glioma stemness has been linked to a number of transcription 
factors, such as S0X2 (Gangemi et al, 2009), 0L1G2 (Ligon 
et al, 2007), and BMll (Facchino et al, 2010). S0X2 and BMIl 
are targets of EMT activators, in particular of ZEBl (Wellner 
et al, 2009). 

ZEBl is an inducer of EMT, transcriptional repressor of cell- 
adhesion molecules, miRNAs — particularly the miR-200 family 
—and cell polarity- associated genes (Brabletz & Brabletz, 2010; 
Wellner et al, 2009). It has emerged as one of the master 
regulators for metastasis (Brabletz & Brabletz, 2010) and plays a 
critical role in tumour initiation at distant sites (Wellner 
et al, 2009). Edwards et al (2011) showed induction of ZEBl 
through the tumour microenvironment in glioma, and related 
ZEBl expression to repression of E-cadherin and thus invasion. 
Given the low prevalence of E-cadherin in glioma (Utsuki 
et al, 2002) , we asked whether ZEBl could affect glioma invasion 
through other mechanisms. We further aimed to elucidate if 
ZEBl regulates other, typically EMT-related processes in brain 
cancer, such as potential for recurrence and therapy resistance. 
Expression of ZEBl has been observed in chemoresistant cells in 
cancers outside the CNS (Li et al, 2012; Wang et al, 2009), but 
whether this relation is correlative or causative is thus far 
unknown. 

Our study provides a systematic analysis of the functions of 
ZEBl in glioblastoma pathobiology. Specifically, we address 



how ZEBl exerts its effects on key malignant processes in 
glioma, i.e. invasion, tumourigenesis and therapy resistance. 

RESULTS 

ZEBl is expressed in invasive glioblastoma cells 

Using three cell lines generated from primary glioblastoma 
specimens (hGBM LO, LI and L2 (DeleyroUe et al, 2011; Piccirillo 
et al, 2006)), we observed varying degrees of tumour invasion 
in xenograft models (Fig lA). In all cases, invasion followed the 
brain architecture, with tumour cells migrating along the 
subcortical white matter across the midline (Fig 1, arrowheads). 
We hypothesised that factors associated with EMT might also 
govern brain tumour invasion. Therefore, we studied the 
expression of EMT-associated factors in our cell lines, and found 
that ZEBl, ZEB2, Twist 1 and Engrailed 1, as well as N-cadherin are 
expressed across all three lines (Fig IB). In contrast. Snail, Slug and 
E-cadherin were found in only one line (Supporting Information 
Fig SI A) . Analysis of these factors in the TCGA dataset (The Cancer 
Genome Atlas Research Network, 2008) revealed that only ZEBl 
showed a significant correlation with patient survival (Supporting 
Information Fig SIB). Hence, we chose to focus on this protein. 

ZEBl is preferentially localised at the invasion front in tumour 
xenografts (Fig IC, Supporting Information Fig SIC). We observed 
that N-cadherin expression is confined to the tumour mass, 
tapering off towards the invasion front, and is absent in invasive 
cells (Fig ID) . Notably, an inverse relationship between N-cadherin 
expression and invasion has been described (Asano et al, 2004), 
and the importance of cell-cell interactions for EMT processes and 
plasticity has recently been discussed (Thompson & Haviv, 2011). 
Beta-catenin is known to induce expression of ZEBl (Kahlert 
et al, 2012; Schmalhofer et al, 2009), but we did not observe 
nuclear accumulation of beta-catenin in our samples. Instead, 
immunoreactivity was confined to cell membranes within the 
tumour core, and absent at the invasion front, corroborating 
reduced cell-cell contacts during invasion (Supporting Information 
Fig SID). Because the ZEBl-positive population was coincident 
with the distal-most, invading tumour cells, we hypothesised that 
ZEBl may have a regulatory role in glioma invasion, and attempted 
to determine its specific functions in glioblastoma. 

ZEBl knockdown reduces invasion and chemoresistance 

To address its functions in invasion we knocked down ZEBl 
expression in glioblastoma cells. Three short hairpin RNA 
constructs (including one previously described (Wellner 
et al, 2009)) were used to knockdown ZEBl in glioblastoma cells 
(Fig 2A) . This resulted in increased localisation of N-cadherin and 
beta-catenin at the membrane (Supporting Information Fig S2A), 
as well as in higher proliferation (Supporting Information Fig S2B). 
Xenograft tumours of ZEBl knockdown cells (shZEBl) were non- 
invasive, and consisted only of an expansive tumour mass 
(Fig 2B) . Of note, invasion correlated with ZEBl expression levels, 
and was highest in the LI and lowest in the LO primary line (Fig 2B, 
Fig lA and B). Immunofluorescence demonstrated an inverse 
relationship between ZEBl and N-cadherin at the invasion front 
of control tumours, while membrane-associated expression of 
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Figure 1. ZEBl is expressed in invasive tumour cells. Source data is available for this figure in the Supporting Information. 

A. Primary glioblastoma (hCBM) cell lines LO, LI and L2 show varying degrees of tumour mass (asterisks) and invasion (arrowheads). Images are inverted black 
and white micrographs of fluorescent immunostainings for human-specific nestin (n = 5 animals each). 

B. EMT factors ZEBl, ZEB2, Twistl and En-1, as well as N-cadherin are expressed in invasive hGBM cell lines. Additional markers are shown in Supporting 
Information Fig lA. 

C. ZEBl is highly expressed at the invasion front of tumours, but not within the tumuor mass (dotted line separates areas of tumour mass and invasion; single 
image projection from confocal z-stack, scale bars 20 ixm). Bar graphs depict ratios of mean fluorescence intensity at randomly selected tumour core and edge 
areas (n = 5 each, Mann-Whitney L/-test). 

D. While the tumour mass is immunopositive for N-cadherin, invasive cells are negative. Scale bars 20 ixm. 



N-cadherin was retained up to the sharply defined borders of the 
shZEBl -derived tumour mass (Fig 2C). 

Since EMT has been associated with therapy resistance 
(Singh & Settleman, 2010), we next tested whether ZEBl affects 



response of tumour cells to the standard of care drug TMZ. A cell 
viability analysis revealed that ZEBl knockdown significantly 
increased TMZ sensitivity in vitro (Fig 2D) . To validate whether 
ZEBl promotes chemoresistance in vivo, we treated tumour- 
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Figure 2. ZEBl knockdown reduces invasion and chemoresistance. Source data is available for this figure in the Supporting Information. 

A. Short hairpin constructs against ZEBl prominently reduce its protein levels. shZEBl targets ZEBl protein in all three primary cell lines. 

B. ZEBl knockdown results in reduced tumour invasion in all cell lines (n = 5 animals per group, Mann-Whitney L/-test). Threshold images used to determine 
invasion index are shown in Supporting Information Fig S2C. 

C. Control tumours show expression of ZEBl, but not N-cadherin, at the invasion front, while tumours derived from ZEBl knockdown cells are characterised by 
high levels of N-cadherin and no invasion. Scale bar 20 ixm (applies to all panels). 

D. ZEBl knockdown results in significantly greater sensitivity to TMZ in MTT cell viability assays (n = 8, one-way ANOVA) and in xenografted tumours (n = 5 
animals each, log-rank test, arrow indicates TMZ treatment). Increased ZEBl expression at the edges of control tumours is associated with expression of 
MGMT. ZEBl knockdown tumours do not stain for MCMT. Scale bars 50 ixm in left and 10 ixm in other panels. 
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bearing animals grafted with either shZEBl or shGFP cells with 
one cycle (five doses) of clinically relevant TMZ concentrations 
(20 mg/kg (Zhou et al, 2007)). This resulted in a significantly 
increased survival time of shZEBl over shGFP animals (Fig 2D). 
Further, we observed that expression of MGMT, a major 
chemoresistance enzyme (Bocangel et al, 2002), increased 
towards the tumour edges in control tumours, and that ZEBl and 
MGMT co-localise (Fig 2D). MGMT was localised to the 
cytoplasm, which has been described (Ishibashi et al, 1994). 

Together, these data suggest that invasive cells are capable of 
evading chemotherapy and support a function of ZEBl in both 
glioblastoma invasion and chemoresistance. 

Glioblastoma invasion is mediated through ZEBl and ROBO 

Based on our observation of reduced membrane-associated N- 
cadherin in invasive glioblastoma cells, we speculated that 
reduced cell-cell contacts contribute to greater invasiveness of 
ZEBl -positive cells. Curiously, although immunofluorescent 
analysis of xenografted tumours indicated that invasive tumours 
lacked N-cadherin (Fig 2C), Western analysis revealed that 
overall N-cadherin protein levels did not change between ZEBl 
knockdown and controls (Fig 3 A) . Corroborating our immuno- 
blot observations, there was no significant difference in the 
mean fluorescence intensity of N-cadherin between ZEBl 
knockdown and control cells in vitro. Flow cytometry analysis 
of live cells demonstrated that fluorescence intensity as well as 
the number of labelled cells was similar between ZEBl 
knockdown and controls (Fig 3B). Since the antibody used 
was directed against the extracellular epitope of N-cadherin, this 
indicated that N-cadherin was not internalised in invasive cells. 

We therefore hypothesised that the distribution of N-cadherin 
might provide the critical difference for tumour cell adhesion. 
Closer scrutiny of N-cadherin immunostaining confirmed that 
membrane distribution of this protein differed between control 
and ZEBl knockdown cells. In control cells, N-cadherin protein 
stippled the cell membrane evenly, whereas ZEBl knockdown 
resulted in the concentration of N-cadherin to the juxtaposed 
membranes between adjacent cells (Fig 3B). Re-distribution of 
N-cadherin from the area of contact between cells may account 
for greater cell motility and may be achieved if N-cadherin is 
disconnected from its intracellular anchor to the cytoskeleton. 
The axon guidance molecule ROBOl can sever the anchorage of 
N-cadherin to the cytoskeleton (Rhee et al, 2002), thereby 
increasing cellular motility, and is expressed in malignant glioma 
(Mertsch et al, 2008). 

A database search for microRNA binding sites revealed ROBOl 
as a potential target of the ZEBl regulatory loop (Ghosh, 2000; 
John et al, 2004). Therefore, we scrutinised ROBOl as intermedi- 
ary regulator of invasion. Consistently, we found that ZEBl 
knockdown reduced expression levels of ROBOl, while over- 
expression of ZEBl had the opposite effect (Fig 3 A) . Additionally, 
ROBOl is located at the outer membrane of control, but not ZEBl 
knockdown cells (Supporting Information Fig S3 A). Overexpress- 
ing or antagonising miR-200c decreased or increased ROBOl 
expression, respectively (Fig 3 A) . Similarly, these manipulations 
restored normal ROBOl expression levels in ZEBl overexpressing 
and knockdown cells (Supporting Information Fig S3B). Antag- 



onising miR-200c also increased cell migration in vitro (Supporting 
Information Fig S3C). 

Forced expression of ZEBl resulted in greater invasiveness of 
xenograft tumours, with tumour cells covering large distances 
along white matter tracts (Fig 3C, arrowheads). Within tumours, 
ROBOl expression increases towards the invasion front, and is 
inversely proportional to the expression of N-cadherin (Fig 3D). 
We next tested whether interference with ROBOl has direct 
effects on tumour cell migration and invasion. Using three 
inducible shRNA constructs against its sequence, we observed a 
prominent, doxycycline-dependent reduction in ROBOl protein 
expression in two cell lines (Fig 3E). In an in vitro scratch assay, 
doxycycline-induced knockdown of ROBOl prominently inhibited 
cell migration (Supporting Information Fig S3D). Expression of 
ROBO could be rescued by a non-targeted construct (Supporting 
Information Fig S3E). Importantly, shROBOl cells gave rise to less 
invasive tumours in animals that were treated with doxycycline 
(Fig 3E, Supporting Information Fig S3F). Finally, the migratory 
phenotype of ZEBl knockdown cells could be rescued by 
overexpressing ROBOl (Fig 3F, Supporting Information Fig 
S3G), while blocking ROBOl prominently reduced migration of 
ZEBl overexpressing cells (Fig 3G, Supporting Information Fig 
S3G). These data support the notion that ROBOl is regulated by 
ZEBl, and that ROBOl is likewise another potential candidate 
molecule for regulating glioblastoma invasion. 

ZEBl regulates MGMT via miR-200c and c-MYB to promote 
chemoresistance 

We postulated above that EMT-associated factors might govern 
increased chemoresistance of invasive cells. Since ZEBl 
knockdown cells are indeed more sensitive to TMZ (Fig 2D), 
we tried to resolve the underlying mechanism. We speculated 
that ZEBl mediates chemoresistance of invasive cells through 
transcriptional regulation of MGMT, and confirmed reduced 
MGMT levels after ZEBl knockdown in immunoblots (Fig 4A). 
MGMT is unlikely a direct target of the ZEBl loop based on 
sequence analysis for binding sites of microRNAs (Ghosh, 2000; 
John et al, 2004) . However, bioinformatics analysis of the MGMT 
promoter from position -1500 to +10 relative to the transcrip- 
tional start site (TSS; Ghosh, 2000) revealed five potential 
binding sites for the proto-oncogene c-MYB, which is regulated 
by miR-200 (Cesi et al, 201 1) . We found c-MYB depleted in ZEBl 
knockdown cells (Fig 4A), indicating it as a potential 
intermediary between ZEBl/miR-200 and MGMT. 

Targeted expression of ZEBl induced c-MYB and MGMT (Fig 4B, 
Supporting Information Fig S4), and increased TMZ resistance 
in vitro (Fig 4B) . Increasing or decreasing levels of miR-200c elicited 
similar changes in c-MYB and MGMT, as well as chemoresistance 
(Fig 4C). Antagonising miR-200c restored expression of c-MYB 
and MGMT as well as chemoresistance in ZEBl knockdown cells 
(Supporting Information Fig S4A), while overexpressing miR-200c 
depleted expression of c-MYB and MGMT and reduced chemo- 
resistance in ZEBl overexpressing cells (Supporting Information 
Fig S4B). These findings substantiated that ZEBl may mediate 
glioma chemoresistance through miR-200c and c-MYB. 

To further characterise this pathway (ZEBl/miR-200/c-MYB), 
we performed additional experiments to test the role of miR-200c 
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and c-MYB in chemoresistance. Using chromatin immunopre- 
cipitation, we confirmed that c-MYB binds to the MGMT 
promoter in glioblastoma cells (Fig 4D). c-MYB knockdown 
reduced expression of MGMT, which was accompanied by 
reduced TMZ resistance in vitro and in vivo (Fig 4E, Supporting 
Information Fig S4C). In contrast, overexpression of c-MYB 
induced MGMT, and increased chemoresistance. Further, 
chemoresistance and MGMT expression of ZEBl -knockdown 
cells could be rescued by over-expressing c-MYB in these cells 
(Fig 4F, Supporting Information Fig S4D) . MGMT was expressed 
independently of ZEBl in these tumours, supporting c-MYB as 
intermediate regulator of MGMT. In contrast, knockdown of c- 
MYB reduced MGMT expression and chemoresistance of ZEBl 
overexpressing cells (Supporting Information Fig S4E) . Finally, 
chemoresistance of shZEBl cells could be rescued by over- 
expression of MGMT, while knocking down MGMT in ZEBl- 
overexpressing cells increased their chemosensitivity (Fig 4G) . 
Knockdown of ZEBl, c-MYB and MGMT could be rescued by 
non-targeted constructs (Supporting Information Fig S4F) . 

Differences in MGMT expression and chemoresistance have 
been attributed to methylation of the MGMT promoter (Hegi 
et al, 2005) . To elucidate whether methylation differences affect 
MGMT expression in the context of ZEBl, we performed bisulfite 
genomic sequencing analysis on the region from -552 to +120 
relative to the MGMT TSS, and found no significant differences 
in methylation between shGFP and shZEBl knockdown cells in 
any of our cell lines (Supporting Information Fig S5). These 
findings demonstrate that ZEBl induces chemoresistance through 
transcriptional regulation of MGMT via miR-200c and c-MYB. 

The ZEBl-miR-200 loop regulates sternness and tumour 
initiation 

Since increased invasion (Cheng et al, 2011) and therapy 
resistance (Bao et al, 2006; Chen et al, 2012) have been observed 
in cancer stem cells (Venere et al, 2011), and ZEBl is a known 
regulator of stemness and S0X2 in other solid tissue cancers 
(Wellner et al, 2009), we hypothesised that ZEBl has similar 
functions in glioblastoma. We identified miR-200c binding sites 
in the sequences of the stemness-promoting transcription factors 
S0X2 and 0LIG2, as well as the stem cell surface marker GDI 33/ 
PROMl (Supporting Information Fig S2) (John et al, 2004). 
Therefore, we tested whether the ZEBl-miR-200 feedback loop 



influenced expression of these proteins. Indeed, knockdown of 
ZEBl decreased protein levels of S0X2, 0LIG2 and CD133, while 
ZEBl overexpression had opposite effects (Fig 5 A and B). 
Overexpressing or antagonising miR-200c inversely affected 
expression of S0X2, 0LIG2 and CD133, supporting inhibition of 
these factors by miR-200c. In a functional assay for stemness, 
decreasing or increasing ZEBl expression resulted in lower or 
higher sphere-forming frequency, respectively (Fig 5C). Impor- 
tantly, orthotopic grafts of shZEBl cells showed a significantly 
decreased capacity for tumour initiation compared to control or 
ZEBl -expressing cells (Fig 5D). We conclude that ZEBl is a 
regulator of stemness in glioblastoma. 

In summary, we found that invasive cells can be distinguished 
from a xenografted brain tumour mass by their expression of 
ZEBl (Fig 5E). Generally, tumour mass cells are stationary due to 
increased cell-cell adhesion, chemosensitive due to lower expres- 
sion of MGMT, and less capable of tumour formation due to their 
lower stemness. In contrast, ROBOl disconnects N-cadherin from 
the cytoskeleton in invasive cells, thus increasing their motility; 
c-MYB induces expression of MGMT, resulting in higher chemo- 
resistance; expression of S0X2 and 0LIG2 result in greater stemness 
and higher capacity for tumour formation. All these processes 
appear to be governed at least in part by ZEBl and miR-200c. 

Protein levels of ZEBl are predictive of clinical outcome 

Having established a molecular model for the actions of ZEBl in 
glioblastoma, we sought to confirm our studies in specimens of 
glioblastoma patients. Our data identify ZEBl as a regulator of 
malignant glioma invasion, which is further corroborated by 
immunohistochemistry in specimens of invasive brain tumours 
of different grades. ZEBl is associated with invasive glioma cells, 
with an increase in ZEBl-immunoreactive cell numbers in 
higher-grade tumours (Fig 6A) . Since expression differences are 
difficult to interpret in immunoperoxidase staining, we per- 
formed immunofluorescence imaging for ZEBl, ROBOl and 
N-cadherin in glioblastoma specimens. While not all glioblasto- 
mas stained for ZEBl, we found increased expression along the 
tumour invasion front (Fig 6B, Supporting Information Fig S6A), 
in accordance with our findings in xenografts (Fig IC) . Similarly, 
ROBOl was found at the invasion front, and its expression was 
mutually exclusive to that of N-cadherin, which was preferen- 
tially at the tumour core (Fig 6C). To further corroborate the 



^ Figure 3. ROBO is a downstream effector of ZEBl that mediates invasion. Source data is available for this figure in the Supporting Information. 

A. Knockdown of ZEBl has no effect on total protein levels of N-cadherin, but robustly depletes ROBOl. Forced expression of ZEBl does not affect N-cadherin, but 
increases ROBOl levels. Likewise, modulation of miR-200 expression affects expression of ROBOl. 

B. N-cadherin immunostaining detects similar levels on control and ZEBl knockdown cells (MFI: mean fluorescence intensity; n = 11, Mann-Whitney L/-test; 
additional images shown in Supporting Information Fig S3). Flow cytometry shows no difference in N-cadherin positive cell numbers in knockdown versus 
control (extracellular epitope; numbers reflect % positive cells). N-cadherin signal is speckled in control cells, but concentrated at cell-cell contact membranes 
in ZEBl knockdown. 

C. Targeted expression of ZEBl increases tumour invasion (n = 5 animals each, one-way ANOVA; arrowheads point to maximum invasion distances). 

D. ROBOl is expressed at the tumour edge, and is inversely proportional to the expression of N-Cadherin. High magnification shows ROBOl-positive tumour cells 
distal to the invasion front. Scale bars, 10 ixm. 

E. Inducible short hairpin constructs against ROBOl can reduce its protein levels in the presence of doxycycline (Dox). shROBO xenografts show reduced invasion 
in the presence of Dox (n = 5 animals each; one-way ANOVA; arrowheads highlight extent of invasion in control conditions; additional data shown in 
Supporting Information Fig S3). 

F. ROBOl rescues migration phenotype in ZEBl knockdown cells (n = 3 for migration assay; additional data shown in Supporting Information Fig S3). 
C. ROBOl knockdown depletes migration of ZEBl overexpressing cells {n = 3, additional data shown in Supporting Information Fig S3). 
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ZEBl pathway in glioblastoma invasion and chemoresistance, 
we performed an immunoblot analysis in glioblastoma samples. 
ZEBl, N-cadherin, ROBOl, MGMT and Engrailedl are frequently 
expressed in glioblastoma (Fig 6D, Supporting Information 
Fig S6B), confirming our experimental findings. Of note, ZEBl 
was expressed in 40-45% of analysed samples. Immunofluores- 
cent staining of patient tumour specimens revealed a 51 ± 9% 
overlap between ZEBl and MGMT (Supporting Information Fig 
S6C), which corroborated the significant correlation between 
ZEBl and MGMT indicated by Western analysis (Fig 6E, 
arrowheads). We also found, in accordance with reduced 
proliferation rates of ZEBl expressing cells in vitro (Supporting 
Information Fig S2) , that ZEBl and the proliferation marker Ki-67 
were mutually exclusive in patient specimens (3.66 ± 0.21%, 
Supporting Information Fig S6C). 

Additionally, we found that protein expression of ZEBl, but 
surprisingly not MGMT, correlated with reduced survival of 
glioblastoma patients, as well as with shorter duration of successful 
TMZ therapy (Fig 6F). We found no significant correlation of 
ZEBl with patient gender, age, or clinical performance (Supporting 
Information Fig S6D) . However, when analysing changes in KPS 
scores over time, a slim but significant correlation with ZEBl 
was observed (Supporting Information Fig S6D). 

Glioblastoma has recently been stratified into different 
molecular subclasses (Phillips et al, 2006; Verhaak 
et al, 2010). Hence, we asked whether ZEBl expression in 
these tumours is associated with any particular subclass. 
Immunoblot analysis of common pathway aberrations in these 
subclasses (Brennan et al, 2009) and other subclass markers 
(Phillips et al, 2006) revealed clustering of ZEBl in the 
proliferative subclass, while ZEBl negative tumours were 
slightly enriched for mesenchymal and proneural subclasses 
(Fig 6E). ZEBl positive tumours showed a significant enrich- 
ment for PCNA and EGFR, while NFl deletions were 
significantly more prevalent in ZEBl negative specimens 
(Fig 6G, Supporting Information Fig S6E). Notably, induction 
of ZEBl through EGF signalling (Ohashi et al, 2010) and 
correlation of EGFR expression and chemoresistance (Murat 
et al, 2008) has been described. Survival analysis was not 
significant for any subclass when stratifying tumour specimens 



into molecular subclasses, but in all cases ZEBl negative 
tumours showed longer survival than ZEBl positive specimens 
(Supporting Information Fig S6H) . Additionally, the proliferative 
subclass demonstrated the shortest overall survival times, and 
the highest number of ZEBl positive samples, further substanti- 
ating the link between ZEBl and outcome. 

Finally, the TCGA dataset (The Cancer Genome Atlas 
Research Network, 2008) supports the prognostic value of the 
ZEBl/miR-200/c-MYB pathway for glioblastoma, as well as 
significant co-occurrence of ZEBl, MGMT and ROBOl in these 
tumours, and decreased levels of EGFR phosphorylation with 
reduced ZEBl expression (Supporting Information Fig S6G 
and H). 

DISCUSSION 

The findings presented here support that inducers of EMT in 
non-CNS cancers can promote single-cell invasion, chemo- 
resistance and tumourigenesis in glioblastoma. While we 
identify a number of EMT-related factors in primary cell lines 
and tumour specimens, ZEBl seems to dominate these 
processes, in up to 50% of patients with glioblastoma. However, 
the apparent heterogeneity of EMT-associated proteins in our 
samples indicates that these molecules may have redundant 
functions, fine-tune downstream pathways and/or serve specific 
purposes in a temporal and spatially regulated fashion. 

Invasion is a common feature of brain malignancies from 
WHO grade II onwards, and we find ZEBl expressed in invasive 
tumours (grade II-IV) , but not in non- invasive neoplasms (grade 
I, data not shown). Thus, expression of ZEBl seems to be 
associated with infiltrating tumour cells across a spectrum of 
glioma. A number of invasion pathways have been identified 
(Dirks, 2001; Hoelzinger et al, 2007; Lefranc et al, 2005), and it is 
likely that several contribute to glioma invasion. Invasive ZEBl- 
positive cells possess a greater degree of freedom for cellular 
locomotion due to the actions of ROBOl, uncoupling N-cadherin 
from the cytoskeleton (Rhee et al, 2002) . This uncoupling allows 
N-cadherin to diffuse freely across the cell membrane, which 
contributes to the ability of invasive cells to break free from 



^ Figure 4. ZEBl regulates MGMT expression and cliemoresistance tlirougli miR-200c and c-MYB. Source data is available for this figure in the Supporting 
Information. 

A. ZEBl knockdown reduces protein levels of c-MYB and MGMT. 

B. Expression of ZEBl increases protein levels of c-MYB and MGMT, and results in higher resistance to TMZ in uitro. 

C. nniR-200c affects protein expression of ZEBl, c-MYB and MGMT, while antagonising miR-200c has opposite effects. Modulation of nniR-200c affects 
chemoresistance to TMZ (n = 8, one-way ANOVA). 

D. Chromatin immunoprecipitation using anti-MYB antibody demonstrates c-MYB binding to the MGMT promoter in primary glioblastoma cells. Indicated are 
putative c-MYB binding sites in the MGMT promoter (Ghosh, 2000). 

E. Knockdown of c-MYB reduces expression of MGMT, while overexpression of c-MYB increases MGMT levels. Both have no influence on ZEBl expression. 
Resistance to TMZ is significantly reduced by c-MYB knockdown in vitro (n = 8, one-way ANOVA). Ablation of c-MYB increases survival of tumour-bearing 
animals, while overexpression of c-MYB negates TMZ effects (n = 5 animals each, log rank test). Experiments with an additional cell line are shown in 
Supporting Information Fig 4. 

F. Chemoresistance in ZEBl-knockdown cells can be rescued by expression of c-MYB (n = 8, one-way ANOVA). Knockdown of ZEBl and simultaneous expression 
of c-MYB abolish differences in survival times after TMZ treatment (n = 5 animals each, log-rank test). Tumours from ZEBl knockdown, c-MYB expressing 
animals are immunoreactive for MGMT, but not ZEBl. Scale bar, 10 ixm (applies to both images) 

C. Chemoresistance in ZEBl knockdown cells can be rescued by expression of MGMT {n = 8, one-way ANOVA), and chemoresistance of ZEBl expressing cells can 
be ablated by knockdown of MGMT (n = 8, one-way ANOVA). 
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Figure 5. ZEBl regulates stemness and tumour formation. Source data is available for this figure in the Supporting Information. 

A. ZEBl knockdown reduces protein levels of S0X2, 0LIC2 and CD133. 

B. Expression of ZEBl increases protein levels of these proteins. Expression of miR-200c reduces, and antagonising miR-200c increases S0X2, 0LIC2 and CD133 
levels. 

C. Tumour sphere formation is decreased by ZEBl knockdown, and increased by its expression {n = 8, one-way ANOVA). 

D. ZEBl knockdown significantly reduced intracerebral tumour formation, whereas expression of ZEBl resulted in increased tumour formation frequency (chi 
square test). 

E. Schematic model depicts interconnection between invasion, chemoresistance and tumour initiation. Invasive cells express high levels of ZEBl, which 
inhibits the expression of miR-200. ROBOl and c-MYB are consequently expressed and affect cell adhesion (by way of ROBOl activity and the disconnection of 
N-cadherin from the cytoskeleton) and chemoresistance (by way of c-MYB and the induction of MGMT). S0X2 and 0LIC2 mediate stemness and tumour 
initiation capacities in these cells. Non-invasive cells in the tumour mass do not express ZEBl; therefore microRNAs inhibit expression of S0X2, 0LIC2, ROBOl 
and c-MYB. N-cadherin remains tethered to the cytoskeleton and inhibits cell motility. MGMT is not expressed and cells are more sensitive to TMZ. 
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neighbouring tumour cells and infiltrate into the surrounding 
tissue. Tumour infiltration may further be facilitated by a variety 
of cellular migratory mechanisms [e.g. chemotaxis, autocrine/ 
paracrine signalling, etc.). 

In solid-tissue cancers outside the CNS, tumour invasion and 
EMT are driven by an E-cadherin to N-cadherin switch (Thiery 
et al, 2009), characterised by the loss of E-cadherin and a 
concomitant increase in N-cadherin. Here, we observe no 
significant levels of E-cadherin in our sample cohort, and also 
find no change in N-cadherin expression. In contrast, our data 
suggests that increased tumour invasion is mediated by a change 
in N-cadherin dynamics, mediated through ROBOl by ZEBl. Of 
note, there are some cases of glioma where E-cadherin is present 
(Lewis-Tuffin et al, 2010). It is tempting to speculate that the 
subset of glioma in which E-cadherin is abundant has a more 
pronounced epithelial character, and might present as gliosar- 
coma. In these tumours, it is conceivable that the 'classic' E-to-N 
switch is the driving force of tissue infiltration. 

Our data strongly support that invasive cells are also more 
resistant to the current standard of care drug TMZ, rendering 
these cells prime candidates for tumour recurrence. An intricate 
regulatory pathway, including miR-200c, c-MYB and MGMT, 
maintains this resistance. To date, MGMT promoter methylation 
is the most reliable prognostic marker for therapy resistance 
(Hegi et al, 2005) . As ZEBl does not influence methylation in our 
primary cell lines, it is possible that ZEBl protein analysis may 
yield prognostic information that complement MGMT methyla- 
tion data. ZEBl, miR-200c and c-MYB may constitute a novel 
pathway for the chemoresistance enzyme MGMT. We observed 
differences in the subcellular localisation of MGMT between 
xenografts (cytoplasmic) and clinical specimens (nuclear). Of 
note, Ishibashi et al. observed that MGMT is present in both 
fractions, and postulated that cytoplasmic MGMT is translocated 
to the nucleus after nuclear MGMT is depleted (Ishibashi 
et al, 1994). It is therefore possible that the nuclear staining 
pattern in patient specimens is due to previous treatment with 
alkylating agents, resulting in nuclear accumulation of MGMT. 
ZEBl was a better predictor of outcome and therapy response at 
the protein level than MGMT. This may indicate that additional 
pathways regulate MGMT, as we found more specimens MGMT 
positive than ZEBl positive. However, the strong correlation 
between ZEBl and MGMT shows that ZEBl -positive cases are 
highly likely to express MGMT, explaining their poor response to 
TMZ. 

In accordance with the cancer stem cell hypothesis, ZEBl is 
linked to the expression of stemness-associated factors and 
tumourigenesis. Our data indicate that critical stem cell 
regulators, such as S0X2 and 0LIG2, are induced by the 
ZEBl-miR-200 feedback loop in glioblastoma. The presence of 
ZEBl -positive cells at the tumour invasion front strongly supports 
an invasive niche that contains cancer stem cells (Cheng 
et al, 2011; Lathia et al, 2011). As the ZEBl-positive population 
is characterised by high motility and increased chemoresistance, 
these cancer stem cells are a candidate population for tumour 
recurrence. Of note, we have previously identified a more 
quiescent population of glioblastoma stem cells (DeleyroUe 
et al, 2011; Piccirillo et al, 2006), and we observed here that 



increased levels of ZEBl are associated with slower proliferation. 
Whether the ZEBl pathway directly reduces proliferation rates of 
glioblastoma stem cells remains to be tested. 

Protein-level analysis shows that about 45% of human 
glioblastomas express ZEBl, and patients with ZEBl -negative 
glioblastomas have survival benefits that are likely related to an 
improved response to TMZ therapy. We cannot exclude that ZEBl 
protein analysis is affected by sample collection site [i.e. core vs. 
edge), but fluorescence immunostaining detected ZEBl only in 
samples that were found ZEBl positive in immunoblots. Sample 
bias may also affect collections in databases [e.g. TCGA), which 
might explain the comparatively low prevalence of ZEBl (and 
other EMT-related factors) in these datasets. Another possibility 
is post- transcriptional regulation of ZEBl expression, e.g. through 
microRNAs, which may result in divergent levels of mRNA and 
protein. Further studies are required to address these issues. 

ZEBl expression in patient tumour samples is significantly 
enriched in the proliferative subclass (Brennan et al, 2009; 
Phillips et al, 2006), which appears in conflict with the lower 
proliferation of ZEBl -expressing cells in patient samples. A 
potential explanation is that ZEBl-positive cancer stem cells 
generate large numbers of rapidly dividing progenies, which in 
turn drive classification. Indeed, it has recently been observed 
that CD133-positive glioblastoma stem cells generate rapidly 
proliferating, but less invasive and less tumourigenic progenies 
(Chen et al, 2010). This study also found significantly lower 
percentages of CD133-positive cancer cells in mesenchymal 
versus proliferative glioblastoma. 

An enrichment of EGFR amplification has been described in 
the proliferative subclass (Huse et al, 201 1) , which we confirmed 
in our cohort. Given the strong correlation with EGFR 
expression, it is conceivable that EGF signalling induces ZEBl 
in these tumours. Since others have found induction of ZEBl 
through beta-catenin or NF-KB signalling (Edwards et al, 2011; 
Kahlert et al, 2012), it is possible that different pathways activate 
ZEBl in different subclasses. The enrichment of EGFR expres- 
sion in ZEBl positive tumours suggests that tumour treatment 
may be more efficacious via EGFR inhibition in combination 
with TMZ treatment. However, clinical trials employing this 
strategy have shown only limited effects to date (Wick 
et al, 2011). Further research is needed to determine whether 
better patient selection or new pharmacological approaches may 
be more successful and whether ZEBl status may be of 
prognostic value in these cases. 

Our findings establish ZEBl as a regulator of invasion and 
chemoresistance in glioblastoma, and a candidate agent for 
tumour recurrence. The multiple ZEBl -associated regulators of 
brain tumour growth and invasion outlined in this study provide 
potential targets for future therapeutic approaches intervening at 
the level of invasion and/or chemoresistance. 

MATERIALS AND METHODS 

Cell culture 

Tumourcell lines were generated (Piccirillo et al, 2006) and maintained 
(Siebzehnrubl et al, 2009) as described. Briefly, 50,000 cells were 



1206 © 2013 The Authors. Published by John Wiley and Sons, Ltd on behalf of EMBO. 



EMBO Mol Med (2013) 5, 1196-1212 



www.embomolmed.org 



Research Article 

Florian A. Siebzehnrubl et ai 



WHO II 



A 

non-neoplastic 

CD • ■ • f 



B Glioblastoma (n=14) ZEB1 Hoechst 



WHO 



WHO IV 





C R0B01 N-Cadherin Hoechst 




160- 

160- 

160- 

20 - 

50 - 
40 - 



Glioblastoma 



ZEB1 

N-Cadherin 

R0B01 

MGMT 

Engrailed-1 

GAPDH 



Correlation of ZEB1 and MGMT 
in Glioblastoma (n=47) 



^ 2.0 - 



p<0.0001 
r=0.6234 




ZEB1 



0 1.0 2.0 

ZEB1/GAPDH (a. u.) 

MGMT Hoechst ZEB1 MGMT Hoechst 



CO 
LJJ 
N 




GBM patient survival by ZEB1 class 



-•-ZEB1+ (n=21) 
-■ ZEB1-(n=29) 




p=0.0137 



20 40 60 

survival (months) 



GBM patient survival by MGMT class 



-.- MGMT+(n=28) 
-.- MGMT- (n=20) 




p=0.9469 



20 40 60 

survival (months) 



E c 10 
N I 5 

I- p 



Chemotherapy by class 
p=0.0206 p=0.7008 



T 



ZEB1 



M 



MGMT 



Molecular subclasses of GBM 
ZEB1 + ZEB1 - 



40— 

110— 
160— 

80— 
40— 

160— 
40— 

50— 
110— 

40— 
60— 



PCNA 

T0P2A 

EGFR 

CD44 
CHI3L1/ 
YKL40 
NF-1 

0LIG2 

DLLS 

N0TCH1 

PDGF-B 

PTEN 



ZEB class for all subjects 

100 
. 80 
I 60 



20 
0 

ZEB1 - 



ZEB1 class by molecular subclass 

25 
. 20 
£ 15 

D 

'O 

1) 10 
5 
0 

ZEB1 - + - + - + 

Mesenchymal Proneural Proliferative 



M 



[1VI]| M I P II R I M |[R]rP]rM]rM]| P I M I M I P |[R] subclass 





Mesenchymal 


Proneural 


Proliferative 


ZEB1 - 


13(72.2%) 


11 (68.8%) 


4 (28.6%) 


ZEB1 + 


5 (27.8 %) 


5(31.2%) 


10(71.4%) 



p=0.027 



Figure 6. 
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seeded per ml of culture medium (N2, Invitrogen, Carlsbad, CA) in the 
presence of mitogens (20 ng/ml each of EGF and FGF2, Sigma, St. Louis, 
MO). Cells were propagated as spheres and passaged using Accutase 
(PAA, Colbe, Germany) every 7 days. For experiments with adherent 
cells, spheres were dissociated and plated in N2 medium supplemented 
with 1% foetal bovine serum (FBS). For scratch assays, 2 x 10^ cells 
were plated per well of a six-well culture plate coated with poly-i- 
ornithine and laminin (15 |JLg/ml) in N2 containing 1% FBS, and grown 
to confluence overnight. Confluent monolayers were scratched with a 
pipette tip, and imaged at the time of the lesion and 24 h later. For 
sphere formation assays, 1,000 cells were plated per well into 96-well 
culture plates in N2 containing EGF and bFGF. Spheres larger than 
50 |jLm were counted 7 days after plating. 

Flow cytometry 

Immunostaining of cancer cells and flow cytometry was performed as 
described (Deleyrolle et al, 2011; Piccirillo et al, 2006). Briefly, live cells 
were dissociated using PBS and 0.5 mM EDTA and subsequently 
incubated with primary antibody or IgG controls in PBS containing 
0.1% BSAfor 1 h on ice. Following 2 washes in PBS, cells were analysed 
on a BD LSR II (BD Biosciences, San Jose, CA). Data was analysed and dot 
plots were generated using Flowjo Ver. 8.8.7 (Tree Star, Ashland, OR). 

Cell viability assay 

The Methyltetrazolium bromide (MTT) assay was used as indicator of 
cell viability and performed as described (Holsken et al, 2006). Briefly, 
10,000 cells were plated per well into 96-well cell culture plates and 
treated 1 h after plating with varying concentrations of TMZ (ranging 
5 ijlM-S mM, Tocris, Ellisville, MO). Concentration-effect curves for 
TMZ treatment were generated by nonlinear regression analysis as 
described (Holsken et al, 2006). Bar graphs are derived from individual 
concentration measurements, compared to the appropriate controls. 

Knock-down experiments 

Plasmids for knockdown of ZEBl and expression of hsa-miR-200c, as 
well as antago-miR-200c and control sequences are as described 
previously (Wellner et al, 2009). Plasmids for knockdown of c-MYB, 
MGMT and ROBOl were obtained from OpenBiosystems (Lafayette, 
CO). The expression plasmid for c-MYB (Clarke et al, 1988) was a kind 
gift of Dr. J.S. Lipsick (Stanford University). Expression plasmids for 



MGMT and ZEBl were obtained from Origene (Rockville, MD). Cancer 
cells were transfected using Lipofectamine LTX (Invitrogen) according 
to the manufacturer's instructions. Transfected cells were selected 
using puromycin orgeneticin (Sigma) before being used for subsequent 
experiments. 

Animal experiments 

Adult female Fox-Chase SCID mice (Charles River, Wilmington, MA) 
were used for in uiuo tumour transplants. All procedures were 
performed according to NIH and institutional guidelines for animal 
care and handling. After animals were deeply anaesthetised using USP 
grade Isoflurane (Halocarbon, North Augusta, SC), an incision was 
made in the scalp, the skull demonstrated and a hole drilled at the 
coordinates Bregma -0.5 mm anterior and -1.5 mm lateral. A 
Hamilton syringe was lowered 2.5 mm into the burr hole, and 1 |jlI of a 
cell suspension was injected over 5 min before the needle was 
retracted. After the incision was closed with surgical staples the animal 
was allowed to recover before being returned to the cage. Animals were 
transplanted with doses ranging from 1,000 to 100,000 cells, and 
tumour-bearing animals were scored regularly for tumour-related 
symptoms. Moribund animals were anaesthetised and transcardially 
perfused with 4% paraformaldehyde in saline, the brains removed, 
postfixed and prepared for histology. For in uiuo TMZ treatment, 
animals received orthotopic grafts of 150,000 cells. One week after 
transplantation, tumour-bearing animals were intraperitoneally 
injected with 20 mg/kg TMZ in saline (final DMSO concentration 
25%). Animals received five injections over 5 days, corresponding to 
one cycle of TMZ treatment. Animals at endpoint (defined as 
bodyweight loss >20% or observation of severe neurological 
symptoms) were perfused and tumour presence was confirmed 
histologically. For in uiuo analysis of inducible ROBOl knockdown, 
animals received orthotopic grafts of 150,000 cells. One group received 
doxycycline injections i.p. (10 mg/kg in saline) at day 4 post surgery 
and then every other day. 

Patient sample collection 

Tissue specimens from glioblastoma patients were obtained from the 
Florida Center for Brain Tumour Research (FCBTR) and the UF and 
Shands Department of Pathology with approval from the UF 
Institutional Review Board (IRB). Informed consent was obtained from 



^ Figure 6. ZEBl expression is associated witli unfavourable outcome in glioblastoma patients. Source data is available for this figure in the Supporting Information. 

A. ZEBl is weakly reactive in some neurons (arrow) and strongly reactive in scattered glial cells (arrowhead) in non-neoplastic parenchyma, and labels glioma 
cells with increasing frequency and intensity with progressing malignancy (grade II astrocytoma, grade III anaplastic astrocytoma and grade IV glioblastoma; 
inset in grade IV depicts pseudopalisading necrosis; scale bars 10 ixm). 

B. In glioblastoma, ZEBl expression is highest at the invasive edge, which is mirrored by expression of ROBOl. 

C. N-cadherin is most prominent at the core, and virtually absent at the edge (representative images from 14 cases, scale bars 20 ixm). Representative images for 
ZEBl negative tumours are shown in Supporting Information Fig 55. 

D. EMT markers ZEBl and En-1 are expressed in glioblastoma samples, as well as N-cadherin and MGMT. Additional immunoblots for EMT-associated proteins 
are shown in Supporting Information Fig S6. 

E. Expression of MGMT and ZEBl are correlated {n = 47, Spearman Correlation). Within ZEBl-expressing tumours, MGMT is frequently co-expressed in ZEB1+ 
cells (arrowheads; 51 ± 9% co-localisation; n = 14; see also Supporting Information Fig 56). Nuclei are counterstained with Hoechst. Scale bar 10 ixm. 

F. Immunoblotting analyses in 50 patient tumour specimens revealed that approx. 40-45% of glioblastomas express ZEBl. This is correlated with clinical 
outcome (median survival ZEB1+ 10.2 months, ZEBl- 15.7 months, log-rank test). While comparable specimen numbers expressed MGMT, no significant 
correlation with outcome was found. Response to TMZ treatment (average TMZ treatment duration ZEB1+ 2.0 months, ZEBl- 7.7 months, Mann-Whitney 
test) was also significantly extended for ZEB1+, but not for MGMT+ specimens. 

C. Stratification of glioblastoma specimens into molecular subclasses revealed a significantly higher portion of ZEB1+ specimens in the proliferative subclass, 
while more ZEBl- specimens were grouped into the mesenchymal and proneural subclasses. Analysis of individual markers across specimens is presented in 
Supporting Information Fig 56. 
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PROBLEM: 

Glioblastoma has a very poor prognosis, mainly due to its 
infiltrating and therapy resistant cells, which eventually lead to 
tumor recurrence. Outside the CNS, epithelial-mesenchymal 
transition has been identified as an inducer of invasion, 
chemoresistance and metastasis in solid tissue cancers, but the 
extent to which similar processes affect brain cancer progression 
is currently unknown. 

RESULTS: 

We find that in glioblastoma the processes of tumor formation, 
invasion and chemoresistance are not separate entities, but 
intertwined. The EMT activator ZEBl is a master regulator of all 
three elements, indicating that EMT-associated factors are 
involved in brain tumor progression. Additional regulatory 
molecules mediate the individual parts of this triad. Tumor 
formation and stemness are associated with increased expres- 



1 



sion of S0X2 and 0LIG2; ROBOl elicits redistribution of N 
cadherin, which reduces cell-cell adhesion and affords greater 
migratory capacity; c-MYB increases expression of the chemo- 
resistance enzyme MGMT. Thus, ZEBl-positive tumor cells evade 
conventional anti-cancer therapy and are a potential source for 
tumor recurrence. Importantly, glioblastoma patients with 
detectable levels of ZEBl have a poorer prognosis and respond 
less to Temozolomide treatment than patients where ZEBl is 
undetectable. 

IMPACT: 

Our findings indicate that (i) a common pathway drives 
tumorigenesis, tissue invasion, and chemoresistance in glio- 
blastoma; and (ii) ZEBl may have prognostic value in 
glioblastoma patients. 



all subjects and all experiments conformed to the principles set out in 
the WMA Declaration of Helsinki and the NIH Belmont Report. 

Immunohistochemistry and Immunocytochemistry 

Immunostainings were performed as described (Siebzehnrubl 
et al, 2009; Zheng et al, 2006). Paraffin-embedded patient material 
was deparaffinised, followed by heat-mediated antigen retrieval in a 
10 mM citric acid buffer (pH 6.0) or Trilogy (Cell Marque, Rocklin, CA), 
and tissue was subsequently immunostained with 3,3'-diaminoben- 
zidine (Vector Elite Kit, Vector labs, Burlingame, CA) or fluorescence- 
conjugated antibodies using standard protocols. A table of employed 
antibodies, suppliers and dilutions can be found in the Supporting 
Information. 

Image acquisition and data analysis 

Low-power fluorescent images were taken on a Leica DMLB 
epifluorescence microscope (Bannockburn, IL) equipped with a CCD 
camera (Spot Imaging Solutions, Sterling Heights, Ml). To obtain full 
images of brain sections, multiple grey-scale images were acquired per 
section using Spot Advanced software (Spot Imaging Solutions) and 
merged into a full image and inverted into black-on-white images 
using Photoshop CS4 (Adobe Systems, San Jose, CA). Photomerged 
images were imported into Image] and threshold levels were adjusted 
to distinguish tumour from background. Using the wand tool, all 
outlines of positively stained (black) tumour areas were selected in 
each section and the perimeter (line surrounding the tumour) and area 
of the tumour were measured. The wand tool allows an exact 
distinction between black (tumour) and white (parenchyma) regions; 
hence, the measurement of tumour outline and area is unbiased. The 
ratio of the squared perimeter distance over the area {P^/A) was 
calculated and used to compare invasive properties of different 
tumours. Since P^/A is a dimensionless number, the resulting figure is 
termed 'invasion index'. A higher invasion index is indicative of a more 
dissociated tumour, whereas a lower invasion index represents a more 



spherical tumour. High-power images were taken on an Olympus BX- 
81 DSU spinning disc confocal microscope (Olympus, Center Valley, PA) 
and projection images of z-stacks were generated using Slidebook 
(Olympus) software. For mean fluorescence intensity analysis, two 
visual fields within the tumour core or edge were selected at random 
per animal, and a confocal z-stack through the entire section was 
obtained. For each stack, one plane of section was selected and mean 
gray values for each channel obtained using Image], and the ratio of 
average channel intensities for ZEBl and Hoechst was calculated. 

RNA isolation and quantitative real-time PGR 

Total RNA was isolated from tumour sphere or adherent cultures using 
the RNeasy Mini Kit (Qiagen, Valencia, CA) according to the 
manufacturer's instructions. RNA was quantified on a Nanodrop 
Spectrophotometer (Thermo, Wilmington, DE), and 1 ijug of total RNA 
was used for cDNA synthesis as described (Siebzehnrubl et al, 2009). 
Twenty-five nanograms of cDNA were used for quantitative PGR using 
the SYBR green PGR master mix (Applied Biosystems, Carlsbad, CA) on an 
ABl 7900HT (Applied Biosystems) as previously described (Siebzehnrubl 
et al, 2009). Expression levels of ZEBl were quantified in triplicate 
relative to beta-actin using the AACt method. Primer sequences and 
amplification times are described elsewhere (Wellner et al, 2009). 

Protein isolation and Western blotting 

Proteins were isolated from cancer cell cultures and primary tumour 
specimens as described (Siebzehnrubl et al, 2009). For Western 
blotting, 5-40 [jug of denatured protein were loaded on 4-12% Bis-Tris 
reducing gel (Invitrogen), separated and blotted onto a PVDF 
membrane (iBlot, Invitrogen). Blots were blocked and probed with 
respective primary and secondary antibodies (see Supporting Informa- 
tion) as described (Siebzehnrubl et al, 2009), and developed using the 
ECL Plus kit (Amersham, Piscataway, NJ) on a FluorChemQ Multi Image 
III (Cell Biosciences, Santa Clara, CA) and Alphalnnotech software 
version 1.0.1.1. Band densitometry was performed using Image]. 
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Chromatin immunoprecipitation 

Chromatin was cross-linked and isolated from 2.0 x 10^ cells in 
sphere culture using the Simple ChIP kit (Cell Signaling, Danvers, MA) 
according to the manufacturer's instructions. Following immunopre- 
cipitation with anti-MYB and IgG control antibodies, the MGMT 
promoter region was amplified with specific primers (primer sequences 
and PCR conditions available upon request). 

Bisulfite genomic sequencing 

DNA was prepared and bisulfite converted as previously described 
(Pardo et al, 2011). Bisulfite converted DNA was amplified using 
MGMT-spec\f\c primers (sequences available upon request). PCR 
products were gel-extracted and cloned into pGEM T-easy vector 
(Promega, Madison, Wl) followed by transformation of TOPIO cells 
(Invitrogen). Individual clones were sequenced and data was analysed 
using Sequencher for sequence alignment and MethylViewer to assign 
site-specific methylation status (Pardo et al, 2011). Molecules with 
<97% conversion efficiency were excluded. 

Statistical testing 

Statistical analyses were performed in GraphPad Prizm 5.0 (GraphPad 
Software, La Jolla, CA). Statistical tests are indicated in the text. In all 
analyses, a p-value < 0.05 was deemed significant. We used the 
D'Agosti no-Pearson test for normal distribution of values. A Bonferroni 
Multiple Comparison test was applied to all ANOVA analyses post-test. 
Observers were blinded to the patient data (including survival time) when 
performing tumour sample protein analyses. We used the R statistical 
software package (V. 2.15.0) to calculate descriptive statistics, create 
graphics and perform all analyses of patient specimens. To evaluate the 
possible associations between molecular subclass (mesenchymal, 
proneural or proliferative) and the other variables in the study (age, 
gender, mortality KPS, ZEBl and MGMT), we used chi-square. Fisher's 
exact or Kruskal-Wallis tests, as appropriate. We used the Kaplan-Meier 
method to perform survival analyses comparing groups classified by ZEBl 
or MGMT level, and to create survival plots. To estimate the effects of ZEBl 
and MGMTwhen controlling for age and gender, we used Cox Proportional 
Hazards models. In all survival analyses, the outcome variable was time 
from start of treatment until death. Subjects still alive at the time or 
analysis and subjects lost to follow-up were considered censored. 
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